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ABSTRACT 

Utilizing a one-dimensional transport theory for heavy ion propagation, 
evaluations of typical energy-independent transport coefficient approximations 
are made by comparing theoretical depth-dose predictions to published experi- 
mental values for incident 670 MeV/necleon Ne beams in water. Results are 
presented for cases where the input nuclear absorption cross sections, or in- 
put fragmentation parameters, or both, are fixed. The lack of fragment charge 
and mass conservation resulting from the use of Silberberg-Tsao fragmentation 
parameters continues to be the main source of disagreement between theory and 
experiment. 
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INTRODUCTION 

As the era of the Space Transportation System progresses toward the 
development of a permanently manned Space Station, possible lunar bases and 
manned Mars missions, long-term exposure of astronauts and spacecraft equip- 
ment to large fluences of galactic heavy ions Indicates a need to investigate 
methods of shielding from these high-energy ions. Because of significant, 
technological barriers, active radiation shielding methods involving electro- 
magnetic fields do not appear to offer viable alternatives to bulk, passive 

shielding methods in the near future (To 83). To properly evaluate passive 
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shield requirements, a comprehensive theory describing the interaction and 
propagation of galactic heavy ions, and their subsequent reaction products , in 
the spacecraft structure and inhabitants is required. In previous work. 

(Wi 84), an energy-dependent transport theory for heavy-ion beams was pre- 
sented and comparisons with an experimental Bragg average ionization curve 
were made. The agreement between theory and experiment was excellent when the 
Si lbe rberg-Tsao fragmentation parameters (Si 77) were renormalized to conserve 
fragment mass and charge* Because of excessive computer resource require- 
ments, however, a fully energy-dependent transport calculation for the com- 
plete galactic heavy ion radiation spectrum is impractical. Typically, 
galactic cosmic ray shielding calculations (Le 83; Wi 86a) utilize energy- 
independent approximations for the input transport coefficients in order to 
reduce the required computer resources. In the work described herein, typical 
energy-independent approximations are used in the fully energy-dependent 
transport theory and the resultant depth-dose predictions compared with pub- 
lished experimental values for incident 670 MeV/nucleon Ne beams in water 
in order to assess the levels of inaccuracy introduced by these input 
approximations . 


DEPTH-DOSE EXPRESSIONS 

From before (Wi 84), the transport equation in the straigh-ahead approx- 
imation and neglecting target secondary fragments is written as 
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where $j(x,E) is the flux of ions of type j with atomic mass Aj at x 
moving along the x-axis at energy E in units of MeV/nucleon, Oj(E) is the 
corresponding macroscopic nuclear absorption cross section, ?j(E) is the 


specific stopping power (change in E per unit distance), and ®j k (E) is the 
fragmentation parameter of ion j produced in collision by ion k. The 
transport equation (1) is solved by the method of characteristics using an 
iterative procedure. The resultant series is 

am 

♦ (x,E) - l * (a) (x,E) (2) 

J a»o ^ 

where a- I, 2, 3, denotes the various generations of reaction products, 

and a - 0 is the primary beam. Equation (2) can be used to evaluate the 
dose as 


D(x) - l / dE S . (E) * (x,E) (3) 

j 0 J J 

where a water density of 1 g/cm^ is assumed. Expressions for the doses due to 

primary, secondary, and tertiary ions are given in equations (9), (12) and 

(15) of our previous work (Wi 84) and are not repeated here. The previous 

calculations used fully energy-dependent values for the transport coefficients 

Sj(E), 0j(E), and Mj^(E). The experimental and calculated doses for 
20 

670 MeV/necleon Ne beams in water are reproduced in Fig. 1. The curve 
labelled ST denotes the calculation using the Silberberg-Tsao fragmentation 
parameters (Si 77) augmented by the light fragment production cross sections 
of Bertini (Be 68). The curve labelled VR denotes the calculation using the 
renormalized fragmentation parameters which conserve fragment mass and charge. 
The multiplicative renormalization factor is 

F - [(Z A )/(Z A ) 1 1/2 
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where Z g and A g are the total fragment charge and mass obtained from the 
semiempirical formulae (Si 77), and Z and Ap are the incident projectile 
ion charge and mass* Typical differences between theory and experiment before 
the Bragg peak are 1-3% for the VR calculations and up to 20% for the ST cal- 
culations. Beyond the Bragg peak the differences are ~25% for the VR calcula- 
tion and nearly a factor of 3 for the ST calculation. 

INPUT PARAMETER APPROXIMATION RESULTS 
Nuclear Absorption Cross Sections 

Typical galactic cosmic ray transport calculations utilize energy- 
independent absorption cross sections, a i j » obtained from some form of the 
Bradt-Peters parameterization (Br 50) 



where r Q and 6 are energy-independent parameters which have been fitted to 
a particular set of cross section data and A^ and Aj are the mass numbers 
of the colliding nuclei. While certainly adequate for high energies where the 
cross sections are nearly asymptotic, significant differences exist (To 86), 
at energies below 2 GeV/nucleon, between experimental data/detailed theoreti- 
cal formalisms and the values predicted by equation (5). 

To test the sensitivity of the dose predictions to the absorption cross 
section energy dependence, the were fixed at their 2 GeV/nucleon values, 

which are representative of the asymptotic results obtained from equa- 
tion (5). The input fragmentation parameters used in the calculations were 
the fully energy-dependent ones. The results are displayed in Fig. 2 as the 
ratios of calculated to experimental doses. For the renormalized fragmenta- 
tion parameter predictions (label VR) the calculated dose is underestimated by 


< 10% before the Bragg peak and by up to 35% beyond the Bragg peak. For the 
unrenormalized fragmentation parameters (label ST) the calculated dose is 
underestimated by up to 33% before the Bragg peak and by almost a factor of 4 
beyond the Bragg peak. 


Nuclear Fragmentation Parameters 

Aside from the use of energy-independent absorption cross sections, 

another possible simplification to the heavy ion transport problem is the use 

of energy-independent fragmentation parameters. To test this approximation, 

20 

dose calculations for the Ne beam in water were performed using fragmenta- 
tion parameters fixed at the values applicable to the incident beam 

energy of 670 MeV/nucleon. The absorption cross sections were fully energy- 
dependent. The results are displayed in Fig. 3 as the ratios of calculated to 
experimental doses. For the VR fragmentation parameters, the calculated dose 
is within 3Z of the experimental dose in the region before the Bragg peak and 
generally within 10Z beyond the Bragg peak. For the ST fragmentation param- 
eters, the calculated dose underestimates the experimental dose by up to 20Z 
before the Bragg peak and by a factor of two beyond it. Thus, as long as 
fragment charge and mass are conserved, the use of energy-independent 
fragmentation parameters may be reasonable. Recently, an energy-independent 
fragmentation model, which conserves fragment charge and mass without 
renormalization, has been developed (Wi 86b) for use in heavy ion transport 
studies and will be incorporated into the existing transport codes in the near 
future . 


Fragmentation Parameters and Nuclear Absorption Fixed 
Finally, we test the sensitivity of the transport calculations to the 


combined use of energy independent fragmentation parameters and nuclear 



absorption cross sections. The f ragmentation parameters were fixed at their 
values for an incident energy of 670 MeV/nucleon and the nuclear absorption 
cross sections were fixed at the asymptotic values for 2 GeV/nucleon incident 
beam energy. The results for 670 MeV/nucleon ^ u Ne in water are displayed in 
Fig. 4 as the ratio of calculated to experimental dose as a function of ab- 
sorber depth. In the region before the Bragg peak, the VR calculation under- 
estimates the dose by up to 13% whereas the ST calculation underestimates by 
up to 35%. Beyond the Bragg peak, the doses are underestimated by < 28% for 
the VR calculation and by factors of 2-3 for the ST calculation. 

Stopping Power 

For all calculations performed in this parameter study, the fully energy- 
dependent stopping powers obtained by summing the electronic and nuclear con- 
tributions according to the methods described elsewhere (Wi 84) are used. 
Results obtained with a simplified power law range-energy relation are pre- 
sented elsewhere (Wi 86a) and will not be discussed in this work. 

CONCLUSIONS 

Utilizing energy-independent transport coefficient approximation, depth- 

20 

dose calculations for 670 MeV/nucleon Ne beams in thick water absorbers have 
been performed and comparisons to experimental data made. The results indi- 
cate that the accuracies of the predicted doses are least affected by the use 
of the energy-independent fragmentation parameters, if they are renormalized 
to conserve charge and mass. The typical error is less than 10% If fully 
energy-dependent absorption cross sections are used. The use of energy- 
independent absorption cross sections resulted in dose underestimates of up to 
35% for the calculations utilizing renormalized fragmentation parameters. 


Clearly, however, the lack of charge and mass conservation in the Silberberg- 
Tsao fragmentation parameters introduces the most significant errors into the 
calculations* Without renormalizing these parameters to conserve fragment 
mass and charge, the predicted doses underestimate the measured doses by 
20-40Z before the Bragg peak and by a factor of 2-4 after the Bragg peak. 
Finally, it should be emphasized for radiation protection studies that care 
must be exercised in interpreting these results. Because of significant 
differences in radiobiological equivalences (RBE) for the various components 
of the mixed radiation field in the absorber, large errors in biological dose 
predictions may exist even though reasonable agreement with the experimental 
Bragg curve is obtained. Hence, it is imperative that transport theories be 
capable of providing accurate predictions of the radiation field composition 
and not just approximate estimates of the total physical dose. 
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Figure Captions 
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Figure 1.- Energy deposition in water by 670 MeV/nucleon Ne ions. Fully 
energy-dependent calculations using Si lberberg-Tsao (ST) and 
renormalized (VR) fragmentation parameters are displayed. Also 
displayed are experimental results from Lawrence Berheley 
Laboratory. The dose relative to zero absorber is displayed on the 
ordinate. 

Figure 2.- Ratio of calculated to experimental doses, as a function of depth 

o n 

in water, for a 670 MeV/nucleon Ne beam. The calculations 
utilized energy-dependent renormalized (VR) and unrenormalized (ST) 
Silberberg-Tsao fragmentation parameters, and energy-independent 
nuclear absorption cross sections. For reference, the Bragg peak 
location is labelled (BP). 

Figure 3.- Ratio of calculated to experimental doses, as a function of depth 

9 O 

in water, for a 670 MeV/nucleon u Ne beam. The calculations used 
energy-independent renormalized (VR) and unrenormalized (ST) 
Silberberg-Tsao fragmentation parameters. The nuclear absorption 
cross sections were fully energy-dependent. For reference, the 
Bragg peak location is labelled (BP). 

Figure 4.- Ratio of calculated to experimented doses, as a function of depth 
in water, for a 670 MeV/nucleon ^ u Ne beam. The calculations used 
both energy-independent nuclear absorption cross sections, and 
Silberberg-Tsao (ST) and renormalized (VR) fragmentation 
parameters. For reference, the Bragg peak location is labelled 
(BP). 
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